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ABSTRACT 

This paper describes the results of computational studies on a single cylinder four stroke cycle 
spark ignition engine using Hydrogen as an alternative fuel to petrol. The simulation is done in 
the professional engine simulation software from AVL Austria named as BOOST. The modeling 
methodology involves the use of first law of thermodynamics for engine as an open system when 
valves are open and engine as a closed system when valves are closed. To include the effect of 
gas exchange in the intake and exhaust manifolds when the valves are open, the modeling is 
done using Navier -Stokes equations for manifolds. The design parameters are fixed by engine 
geometry. A matrix was prepared for the operating variables to carry out the simulation. First the 
data was used as per petrol engine needs, and results for engine performance and emissions 
characteristics were generated. The operation was revised with data for proposed hydrogen 
engine and its performance and emissions characteristics were studied. The software gave 
successful results in both the cases. It was observed that the power output was decreased with 
hydrogen as an alternative fuel due to lower volumetric efficiency, the brake specific fuel 
consumption was comparable in both cases due to higher calorific value of hydrogen. The 
emissions were also reduced with hydrogen as fuel because of no carbon atoms as compared to 
petrol fuel. No CO emission was produced with hydrogen, there was reduction in NOx emission 
and also there was drastic reduction in hydrocarbon emissions with hydrogen as fuel. It is 
proposed that hydrogen can be successfully used in petrol engine as an alternative future fuel. 
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INTRODUCTION 

Computer simulation studies helps to predict the behavior of the engine in the petrol and 
Hydrogen fuel modes. We prepare the models for Petrol and Hydrogenfuel modes for the engine 
systems and feed the actual data corresponding to the design and operating conditions of the 
system. It helps to simulate the results without actually performing experiments. Thus a lot of 
money and time is saved. Moreover we can simulate and compute those results which are very 
difficult to be measured experimentally. Favorable computed results pave the way for further 
experimental investigations. 

The objectives are to investigate the feasibility of Hydrogenas alternative fuel in petrol engines. 
It is also intended to choose Hydrogen as fuel because it has no carbon atoms as compared to 
petrol and therefore will produce minimal pollution. Table 1 at the end gives physico-chemical 
properties of Hydrogen and petrol which help us to investigate the feasibility of using Hydrogen 
as an alternative fuel to petrol. 

THEORETICAL BASIS 

The theoretical background including the basic equations for all elements used in the present 
model is summarized below to give a better understanding of the program. 
The Cylinder , High Pressure Cycle, Basic Equation. 

The calculation of the high pressure cycle of an internal combustion engine is based on the first 



law of thermodynamics: 


d(m,c.u) 
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da 
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(Eq.l) 



= change of the internal energy in the cylinder. 
= piston work. 



da 



= fuel heat input. 



y^dQ^ _ ^ eat losses 
da 
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= enthalpy flow due to blow-by 
= blow-by mass flow 

The first law of thermodynamics for high pressure cycle states that the change of internal energy 
in the cylinder is equal to the sum of piston work, fuel heat input, wall heat losses and the 
enthalpy flow due to blow-by. 

Eq.l is valid for engines with internal and external mixture preparation. However the terms 

which take into account the change of gas composition due to combustion, are treated differently 

for internal and external mixture preparation. 

For internal mixture preparation it is assumed that 

The fuel added to the cylinder charge is immediately burnt. 

The combustion products mix instantaneously with the rest of cylinder charge and thus form a 
uniform mixture. 

As a consequence, the Air-Fuel ratio of the charge diminishes continuously from a high value at 
the start of combustion to the final value at the end of combustion. 

In order to solve this equation, models for the combustion process and the wall heat transfer, as 
well as the gas properties as a function of pressure, temperature, and gas composition are 
required. 

Together with the gas equation 

p c = _L .m c .R .T c (Eq.2) 

Establishing the relation between pressure, temperature and density, Eq. 2 for in-cylinder 
temperature can be solved using a Runge-Kutta method. Once the cylinder gas temperature is 
known, the cylinder gas pressure can be obtained from the gas equation. 



Combustion Model 

The following equation for the stoichiometric air requirement specifies how much air is required 
for a complete combustion of 1 kg fuel: 



L st = 137.85 .( 



+ - 



h 



+ - 



o 



12.01 4.032 32.06 32.0 



) [kg Air/kg Fuel] 



-(Eq.3) 
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For lean combustion, the total heat supplied during the cycle can be calculated from the amount 
of fuel in the cylinder and the lower heating value of the fuel. The lower heating value is a fuel 
property and can be calculated from the following formula: 

H u = 34835 x +93870 . h +6280 . n +10465 . s -10800 . o -2440 . w [kj/kg] (Eq.4) 

In rich air fuel mixture combustion, the total heat supplied during the cycle is limited by the 
amount of air in the cylinder. The fuel is totally converted to combustion products even if the 
amount of air available is less than the amount of stoichiometric air. 

Heat Release Approach. 

The vibe function is used to approximate the actual heat release characteristics of an engine: 

d^ = _^_ (m+1) y m e -a.y(m + l) (Eq 5) 

da Aac 

dx = ^ (Eq.6) 

Q 

y = a-- (Eq.7) 

Aac 

The integral of the vibe function gives the fraction of the fuel mass which was burned since the 
start of combustion: 

x= [(— .da) = l- e " a - y(m+1) (Eq.8) 

J da 



Gas Exchange Process , Basic Equation 

The equation for the simulation of the gas exchange process is also the first law of 
thermodynamic s : 



d{m,c.u) 



p c .dV ^ dQw 



drm 



dm- 



-(Eq.9) 



da da da da.fa da.he 

The variation of the mass in the cylinder can be calculated from the sum of the in-flowing and 

out-flowing masses: 

drric 
da 



^ dm ^ drru 
da da 



-(Eq.10) 



Piston Motion 
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Piston motion applies to both the high pressure cycle and the gas exchange process. 

For a standard crank train the piston motion as a function of the crank angle a can be written as: 



s= (r+l).cos\|/-r.cos(\|/+a)-l. Jl - {y .sin<y + a) - -} 2 (Eq. 11) 



\|/ = arcsin( ) 

r + l 



-(Eq.12) 



Heat Transfer 

The heat transfer to the walls of the combustion chamber, i.e. the cylinder head, the piston, and 
the cylinder liner, is calculated from: 

Q wi = Ai .a w . (T c -T wi ) (Eq.13) 

In the case of the liner wall temperature, the axial temperature variation between the piston TDC 
and BDC position is taken into account: 

1 — cx 

l-e 



T, =T 



L,TDC 



X.C 



-(Eq.l4) 



Tl.tdc 
c = ln{- } 



-(Eq.l5) 



Tl, BDC 

For the calculation of the heat transfer coefficient, the Woschnil978 heat transfer model is used. 



Woschni Model 

The woschni model published in 1978 for the high pressure cycle is summarized as follows: 



1 on rk-0.2 0.8 ^ -0.53 

aw = \3Q.D .p c Tc 



Cl.Cm + C2. \p c - pc, o) 



0.8 



-(Eq.16) 



pci.Vci 

CI = 2.28+0.308.cu/cm 
C2 = 0.00324 for DI engines 
C2 = 0.00622 for IDI engines 
For the gas exchange process, the heat transfer coefficient is given by following equation: 

aw = l3Q.D^ 2 .pc *Tc^ 53 .(C3.Cm) 0S (Eq.17) 

C 3 = 6.18+0.417.c u /c m 
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Fuel Injector 

The fuel injector model is based on the calculation algorithm of the flow restriction. This means 
that the air flow rate in the fuel injector depends on the pressure difference across the injector 
and is calculated using the specified flow coefficients. In addition, the amount of fuel specified is 
fed into the air flow. 

In the case of carburetor model, the fuel flow is set to a specified percentage of the 
instantaneous mass flow. 

For the injector model, a measuring point must be specified at the location of the air flow meter. 
In this case the mean air flow at the air flow meter location during the last complete cycle is used 
to determine the amount of fuel. As is the case for continuous fuel injection, the fuelling rate is 
constant over crank angle. 

The fuel is added in gaseous form to the pipe flow. No evaporation is considered. 



Pipe Flow 

The one dimensional gas dynamics in a pipe are described by the continuity equation 
dp 8(p.u) 1 dA 

— = p.u. — . — , (Eq.18) 

dt dx A dx 

the equation for the conservation of momentum 

djpMl = _ d(p.u 2 + p) _ p mZ j_ M _ F* ? (Eq l9) 

dt dx A dx V 

and by the energy equation 

d_E = _ d[u.(E + p)] _ u(E + p) X^ + *L . (Eq.20) 

dt dx A dx V 

The wall friction force can be determined from the wall friction factor Xf : 

FR AT II .-^ ! — _| . 

= .p.u.\u\ (Eq.21) 

V 2.D 11 

Using the Reynold's analogy, the wall heat flow in the pipe can be calculated from the friction 
force and the difference between wall temperature and gas temperature: 

— = — — .p\u\.c P .(Tw - T) (Eq.22) 

V 2.D 11 

During the course of numerical integration of the conservation laws defined in the Eq.20, Eq.21 
and Eq.22, special attention should be focused on the control of the time step. In order to achieve 
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a stable solution, the CFL criterion (stability criterion defined by Courant, Friedrichs and Lewy ) 
must be met: 

Ax 



At< 



u + a 



-(Eq.23) 



This means that a certain relation between the time step and the lengths of the cells must be met. 
The time step to cell size relation is determined at the beginning of the calculation on the basis of 
the specified initial conditions in the pipes. However, the CFL criterion is checked every time 
step during the calculation. If the criterion is not met because of significantly changed flow 
conditions in the pipes, the time step is reduced automatically. 

An ENO scheme is used for the solution of the set of non-linear differential equations discussed 
above. The ENO scheme is based on a finite volume approach. This means that the solution at 
the end of the time step is obtained from the value at the beginning of the time step and from the 
fluxes over the cell borders 



RESULTS AND DISCUSSION 



Effect of Speed on Power 

The Fig. 1 below shows the effect of speed on power. It is seen as the speed increases the power 
also increases due to more number of power cycles per unit time. Further it is seen that the power 
developed by petrol as fuel is higher than hydrogen due to higher volumetric efficiency with 
petrol. 



A Monthly Double-Blind Peer Reviewed Refereed Open A ccess Interna tional e-Journal - Included in the International Serial Directories 
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., MJJiBBtffli^ as well as in Cabell's Directories of Publishing Opportunities, U.S.A. 

International Journal of Management, IT and Engineering 

http://www.ijmra.us 




Effect of Speed on Brake Specific fuel Consumption, (bsfc) 

Fig. 2 below shows the effect of speed on brake specific fuel consumption (fuel consumed per 
unit power output ) . It is seen that the operation of hydrogen engine is approximately as 
economical as the petrol fuel.From the physical and chemical properties of petrol and hydrogen 
it is seen that hydrogen has higher heating value but at same time its volumetric efficiency is 
lesser. 
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Fig .2 Effect of Speed on Brake 
specific fuel consumption 
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Effect of Speed on Torque. 

Fig. 3 below shows the effect of speed on torque. It is seen that in each case the maximum torque 
is produced at 2000 r.p.m.. This is because the combustion characteristics and in-cylinder 
pressure development is best at this speed. Petrol produces better torque due to its higher 
volumetric efficiency as compared to hydrogen. 




Petrol 
■*— Hydrogen 



Fig 



2000 4000 6000 8000 
Speed , r.p.m. 

.3 Effect of Speed on Torque 



Effect of Speed on Exhaust Gas Temperature. 

Fig.4 below shows the effect of speed on exhaust gas temperature. It is seen from the graph that 
the exhaust gas temperatures in case of petrol is higher. This is also a clear indication that petrol 
produces higher temperatures due to overall better combustion characteristics. 
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Fig. 4 Effect of Speed on Exhaust Gas 
Temperture 



Effect of Speed on HC Emissions 

Fig. 5 below shows the effect of speed on hydrocarbon emissions. The HC emissions are much 
less with hydrogen fuel. The contribution of HC emissions with hydrogen fuel is mainly due to 
crank case lubrication which circulates the hydrocarbon based lubricating oil into cylinder for 
piston cylinder lubrication. 
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Fig.5 Effect of Speed on HC 
Emissions 
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Effect of Speed on CO Emissions. 

Fig. 6 below shows the effect of speed on carbon monoxide emissions. CO emissions with petrol 
are significant due to many carbon atoms in petrol. No CO emissions are produced with 
hydrogen fuel as there is no carbon atom in hydrogen. The CO emissions with petrol are higher 
at higher speeds as there are more number of power cycles at higher speed. 
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Fig .6 Effect of Speed on CO 
emissions 



Effect of Speed on NOx Emissions 

Fig.7 below shows the effect of speed on NOx emissions. The NOx emissions with petrol fuel 
are higher as higher temperatures are produced with petrol as compared to hydrogen. For the 
formation of NOx emissions presence of oxygen and nitrogen from air at higher temperatures for 
long time is desirable and favourable condition. 
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Fig .7 Effect of Speed on NOx 
Emissions 



CONCLUSIONS 

1. Hydrogen can safely be used in petrol engines as an alternative fuel to conventional petrol. 

2. Pollution formation from engines using Hydrogen as fuel will be much less. 

3. Lesser power is produced by Hydrogen fuel as compared to petrol. 

4. The Hydrogen engine operation is as economical as petrol approximately. 
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APPENDIX-A 
NOMENCLATURE 

a 
A 

A ef f 

Ai 

AFcp 
A a 
c 

c v 



igeo 



CI 

C2 

C2 

Cm 

Cu 

C u 

D 
D 

dm; 

dm e 

d v i 

drriBB 
da 



f 

F R 
h 

h B B 
hi 



speed of sound 
pipe cross-section 
effective flow area 

surface area (cylinder head, piston, liner) 
air fuel ratio of combustion products 
geometrical flow area 
mass fraction of carbon in the fuel 
specific heat at constant volume 
specific heat at constant pressure 
2.28+0.308.cu/cm 
0.00324 for DI engines 
0.00622 for IDI engines 
mean piston speed 
circumferential velocity 
circumferential velocity 
cylinder bore 
pipe diameter 

mass element flowing into the cylinder 
mass element flowing out of the cylinder 
inner valve seat diameter (reference diameter) 

blow-by mass flow 

piston pin offset 

energy content of the gas (=p. cv T + i .p.u 2 ) 

fraction of evaporation heat from the cylinder charge 
wall friction force 

mass fraction of hydrogen in the fuel 
enthalpy of blow-by 
enthalpy of in-flowing mass 
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h e = enthalpy of the mass leaving the cylinder 

H u = lower heating value 

k = ratio of specific heats 

1 = con-rod length 

m = shape factor 

m = mass flow rate 

m c = mass in the cylinder 

m ev = evaporating fuel 

m p i = mass in the plenum 

n = mass fraction of nitrogen in the fuel 

o = mass fraction of oxygen in the fuel 

p = static pressure 

Poi = upstream stagnation pressure 

Pc,o = cylinder pressure of the motored engine[bar] 

Pc,l = pressure in the cylinder at IVC[bar] 

p p i = pressure in the plenum 

p c = cylinder pressure 

P2 = downstream static pressure 

q ev = evaporation heat of the fuel 

q w = wall heat flow 

Q = total fuel heat input 

Qf = fuel energy 

Q w i = wall heat flow (cylinder head, piston, liner) 

r = crank radius 

Ro = gas constant 

s = piston distance from TDC 

t = time 

T = temperature 

Tc,l = temperature in the cylinder at intake valve closing (IVC) 

T c = gas temperature in the cylinder 

T w ; = wall temperature ( cylinder head, piston, liner) 
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T L = 

T L,TDC = 
T L,BDC = 

T w — 

Toi = 

u = 

u = 
V 
V 

VD = 

w = 

X = 

X = 

a = 

a D 
Aa c 
a w 
P 



1 1 nor temperature 

liner temperature at TDC position 

liner temperature at BDC position 

pipe wall temperature 

upstream stagnation temperature 

specific internal energy 

flow velocity 

cylinder volume 

cell volume (A.dx) 

displacement per cylinder 

mass fraction of water in the fuel 

relative stroke (actual piston position related to full stroke) 
coordinate along the pipe axis 
crank angle 



= start of combustion 

= combustion duration 

= heat transfer coefficient 

= density 

= flow coefficient of the port 
crank angle between vertical crank position and piston TDC position 

A,f JT = wall friction coefficient 

At = time step 

Ax = cell length 
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Petrol Engine Specifications 


Bore 


84 mm 


Stroke 


90 mm 


Compression Ratio 


9 


Number of Cylinders 


1 



APPENDIX-C 



Table 1 : Physico-Chemical Properties of Petrol and Hydrogen [2] 



Fuel Property 


Hydrogen 


Petrol 


Formula 


H2 


C4 to C12 


Molecular weight 


2.02 


100-105 


Lower heating value, MJ/Kg 


121 


42.5 


Stoichiometric air-fuel ratio, 


34.3 


14.7 


weight 






Octane number 


130 


80-98 
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